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The screened electron-electron interaction in a multi-band electron system is calcu- 

O 

lated within the random phase approximation and in the tight-binding representation. 



The obtained dielectric matrix contains, beside the usual site-site correlations, also the 
site-bond and bond-bond correlations, and thus includes all physically relevant polariza- 
tion processes. The arguments are given that the bond contributions are negligible in the 
long wavelength limit. We analyse the system with two non-overlapping bands in this 
limit, and show that the corresponding dielectric matrix reduces to a 2 x 2 form. The 
intra-band and inter-band contributions are represented by diagonal matrix elements, 
while the off-diagonal elements contain the mixing between them. The latter is absent 
in insulators but may be finite in conductors. Performing the multipole expansion of the 



bare long-range interaction, we show that this mixing is directly related to the symme- 
try of the atomic orbitals participating in the tight-binding electronic states. In systems 
with forbidden atomic dipolar transitions, the intra-band and inter-band polarizations 
are separated. However, when the dipolar transitions are allowed, the off-diagonal el- 
ements of the dielectric matrix are of the same order as diagonal ones, due to a finite 
monopole-dipole interaction between the intra-band and inter-band charge fluctuations. 
We also calculate the macroscopic dielectric function and obtain an expression which 
interpolates between the well-known limits of one-band conductors and pure insulators. 
In particular, it is shown that the microscopic origin of the so-called self-polarization 
corrections is the on-site interaction which exchanges two electrons at different orbitals, 
combined with a finite tunneling between neighboring sites. 

PACS: 71.10,+x 

Key words: dielectric matrix, random phase approximation, self-polarization cor- 
rections 
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1 Introduction 

The present and the forthcoming |lj work revisit the problem of the dielectric properties 
of multiband metals and insulators. Early calculations || were carried out neglecting the 
discreteness of the lattice in the derivation of the screened Coulomb interaction. The re- 
sult was a dielectric function built adding the intraband and interband polarizabilities [2 
- 4] . The continuous approximation in calculating the Coulomb matrix elements is how- 
ever inconsistent with the multiband assumption, since the latter implies the existence 
of the lattice. Later works attempted to include properly the local electric field effects 
associated to the lattice discreteness, by formulating the problem of dielectric response in 
the representation of the k+G plane waves, where G are vectors of the reciprocal lattice 
|5|, |J. In this representation one comes to the dielectric matrix £(q+G, q+G', uo) which 
is of infinite order for a given q. Since this latter is very difficult to handle in physical 
terms, it was necessary either to make additional approximations [7 - 12] (beside the 
initial random phase approximation (RPA)), or to develop various numerical procedures 
[13 - 15], which in turn usually do not provide a simple physical interpretation of the 
results obtained. 

A particular method in this direction is based on the use of the tight-binding (TB) 
basis in the calculation of the polarizabilities which figure in the matrix elements of 
£:(q+G, q+G',CG>) [8 - 12]. Providing that the atomic orbitals which represent this 
basis are localized enough, one gets to a good approximation a factorizable form of 
dielectric matrix elements, introduced earlier on through the so-called generalized shell 
model by Sinha et al 0. This factorization in principle enabled the analytical inversion 
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of the dielectric matrix, and led to significant improvements in the understanding of 



dielectric properties of simple TB band structures [g, [T|. The whole approach was still 
burdened by the unnecessary parallel use of two basis, i.e. plane waves and TB states. 

The aim of the present paper is to show that the problem of the dielectric response 
is significantly simplified when the TB scheme is used from the outset, in the definition 
of the dielectric matrix itself. The TB approach for the single band case was developed 
before jOj). In the multiband case the elements of the TB dielectric matrix have a 
transparent physical content, i.e. they are defined by the pairs of intraband and/or 
interband transitions and by an index which covers pure site-site correlations as well 
as the site-bond and bond-bond correlations between neighboring lattice points. By 
taking into account these bond correlations we complete our earlier treatment of the 
dielectric response in the TB approach |17 |. Since the atomic orbitals which build the 



TB Bloch functions are well localized (in contrast to usual Wannier functions), it is well 
justified to keep only the bond contributions between first neighbors, so that the index 
mentioned above is of the order of a coordination number for a given lattice. After 
an additional truncation by keeping a finite number of presumably relevant electronic 
bands, the dielectric matrix reduces to a finite, rather low order. 

In order to illustrate the method, we consider here in some detail the simple case 
of two non-overlapping electron bands, allowing for two types of symmetry of the TB 
orbitals which form the empty higher (conducting) band. In all cases the lower band is 
assumed to be partially or fully filled. We calculate all screened TB Coulomb matrix 
elements within RPA, and show that in the long-wavelength limit the problem reduces to 
a 2 x 2 dielectric matrix. The same result was obtained by an alternative analysis via the 



Heisenberg equations for the electronic charge density |L7|. Since the present approach 



gives in addition screened Coulomb matrix elements for all intraband and interband 
processes, it can be conveniently used in the calculation of the screened local fields, such 
as those produced by phonons in ionic multiband metals, and in particular in high T c 
superconductors. 

Since the range of the TB atomic (molecular) orbitals is shorter than that of the two- 



site Coulomb interactions |TJ| , the latter can be expanded into the multipole series, which 
is particulrly well controlled in the long wavelength limit. In this limit one recognizes the 
decisive role of the symmetry of TB orbitals in the multiband dielectric screening. When 
two orbitals have the same parity, the intra-band and inter-band polarization processes 
are decoupled, so that the total dielectric function is the product of the corresponding 
intra-band and inter-band dielectric functions [19 - 21]. A more interesting situation 
occurs when the band symmetry allows for finite dipolar transitions. The finiteness of 
the monopole-dipole coupling in this case leads to qualitatively different results for the 
microscopic dielectric function (defined as the determinant of the dielectric matrix) and 
the TB matrix elements of screened electron-electron interaction. The effects of this 
coupling on the collective modes of the two-band electron liquid will be analysed in the 
paper II [|I|]. 

The multipole expansion also separates the long-range part of the bare Coulomb 
interaction from the on-site contributions. Because of that, we are able to trace the 
microscopic origin of the so-called self-polarization corrections || to the macroscopic 
dielectric function em{% w), which can be easily determined once the TB matrix elements 
of the screened Coulomb interaction are known. 



The paper is organized in the following way. In Sect. 2 we derive the linear system of 
RPA equations for the screened Coulomb interaction and introduce the corresponding 
dielectric matrix. In Sects. 3 to 5 we apply these general results to the two-band model. 
The explicit expressions for the microscopic dielectric function and all matrix elements 
of the screened interaction are derived in Sect. 3. In Sect. 4 we perform the multipole 
expansion and discuss the role of the short-range and long-range contributions to the 
microscopic dielectric function. The macroscopic dielectric function and the origin of the 
self-polarization corrections are considered in Sect. 5. Sect. 6 contains some concluding 
remarks. 

2 The tight-binding formulation of the dielectric matrix 

The Dyson's equation for the screened interaction in the direct space is within RPA 
given by 



V(r,r» = y(r-r')+ j dv, j dx' x V(r - n)n(n, 4 u)V{r' x , r , u), (1) 



where 





is the bubble polarization diagram, and G°(r, r\u) = J2i Gf(r, r , to) is a bare Green's 



function with the contribution from the Z-th band given by 



G?(r,r',u;) = £^k(r)^ k (r') x 



k 




(3) 
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with ipi^r), Ei(k) and n\ being the corresponding Bloch function, band dispersion and 
occupation number respectively. Our aim is to transform the integral equation (|l[) into 
a system of algebraic equations for the TB matrix elements of the screened potential. 
To this end we recall the essential property of the TB basis, namely the weak overlap 
of the atomic orbitals at neighboring crystal sites. Thus, it suffices to keep in various 
matrix elements only the contributions with pairs of orbitals centered at the same and 
neighboring sites, and to neglect all contributions with more distant pairs of orbitals. 
In this respect it is important to compare the TB and the Wannier bases. Although 
the Wannier functions at different sites are exactly orthogonal, they usually have slowly 
converging oscillatory tails, which greatly complicates the evaluation of the one- and 
two-body matrix elements. We therefore continue by strictly using the TB basis. The 
Bloch functions are thus given by 

Here 

S t (k) = £ e^Sffl, (5) 
5^o 



where 



Si(*) = J dv^{v) Vl {v - 8) (6) 

is the direct overlap between nearest neighboring atomic orbitals (pi, and the sum in (H) 
involves only first neighbors. Although the further considerations mainly do not depend 
on the details of the band dispersions Ei(k), we remind that they include, beside <Sj(k), 
also sums of tunneling integrals due to the effective ionic potential. In these sums it 
again suffices to keep only nearest neighbors. 



Let us now multiply eq.(|T]) by e t ^ R+ ' 1 ^^if* 1 (r — K — 81)^*^')^^ (r — K)^^ (r' — 8 2 ), 
where R, 61 and 8 2 denote crystal sites. After the integration with respect to r and r' 
and summation with respect to R one gets the linear system of equations 



EE 

hi' 3 8 3 



x 



x Vi> 3hl3l > 2 (8 3 , 8 2 ] q, u) = V hhfA (<5i, <5 2 ; q) (7) 



-* '-^ x 



for the matrix elements of the screened Coulomb interaction 

V hhl , A {5 u 8 2 ; % u) =Y,e^ R+ ^ Z) Jdr J dv'^r - R - 8 1 )^ 2 (v' 

xy(r,r»^(r-R)^(r'-<J 2 ). (8) 

V, , (61, 62', q) are the corresponding bare matrix elements, when V(t,t',u) in (0) is 
replaced by e 2 /|r — r'|. In this TB representation of two-body interaction one again 
relies on the localized nature of atomic orbitals and retains only those matrix elements 
in which <5i and 62 are either zero or denote nearest neighboring sites. In other words, 
one distinguishes between the site-site (Si = 62 = 0), site-bond (S\ or #2 equal to zero) 
and bond-bond (8\ and 82 different from zero) contributions to the given q-component 
of the Coulomb interaction. In eq.(^) we have also passed to the adiabatic screened 
matrix elements, with the polarization diagrams given by 

n, „ (8 - 8'- a u) = — V e -*MX*-*') 1 x 

lA ' q ' ] N^ 6 [l + 5i(k)][l + 5i,(k + q)] 

niik) -n^jk + q) 

u + E l (k)-E l ,(li + q)+ir } ' { ] 



k 



The coefficient on the left-hand sides of eq. (|7|) , 

^[Suhi'sSs = kiAsi'Ms, -E^isi^^i'^jqJn^ ^a - ^q,^) (10) 

8' 3 



define the dielectric matrix in the TB representation. Note that the sums in eqs.(|7|) and 
fllQf) go over S = and all nearest neighboring sites. The order of the matrix fllQ|) (i.e. 
the number of equations in the system (|7]) for fixed I2, l 2 and £2) is equal to the number of 
pairs of band indices (Mi) multiplied by Z+l, where Z is the number of nearest neighbors 
(i.e. the coordination number for Bravais lattices). Note that the polarization diagram 
(P) vanishes if l-th and l -th bands are both full or empty, so that the number of relevant 
pairs (//') in eqs.((|) and ( |T0D reduces to that representing transitions between (partially) 
full and (partially) empty bands (including the intra-band transtions). Furthermore, the 
order of the matrix ([!(]) reduces to a finite number after a physical truncation of the 
system (^) by which one keeps only contributions from bands which are not too far 
energetically from the Fermi level, and neglects e.g. those from deep core atomic states 
and very high empty bands. 

The determinant of the matrix (|10D is of the central importance for further con- 
siderations, since it carries relevant information about the microscopic dielectric prop- 
erties. Indeed, it enters into the denominators of all screened matrix elements, and, 
furthermore, its zeros define the collective as well as the electron-hole excitations of the 
two-band electron gas [IJ. For these reasons we call it microscopic dielectric function, 
e m (q,u;) = det[e). 

We summarize this Section by comparing the present TB formulation of the dielectric 
matrix with the more usual one, based on the plane wave representation and indexed 
by an infinite number of the vectors of reciprocal lattice. The expressions (|7| - |10[) 



formally look even less convenient, since the dielectric matrix (|H]) is spanned by two 
multiplied infinities, the number of lattice sites (indices 8) and the number of intra- 



and inter-band transitions [indices (it)]- However, as argued above, both these infinities 
are physically reduced to finite orders, which are eventually rather low in standard TB 
systems. No such reduction is possible in the plane wave formulation, so that the order of 
the corresponding dielectric matrix remains physically very large. The reason is that the 
base of reciprocal lattice vectors cannot be simply truncated whenever the polarization 
processes take part on the scale of the unit cell. 



Two-band model 



The further analysis is concentrated on the case of two electron bands which do not 
overlap. To fix ideas we assume that the lower ("valence") band (I = 0) is partially or 
completely full, while the upper ("conducting") band (I = 1) is empty. The correspond- 
ing band dispersions Ei(k) will be specified in the later examples Furthermore, the 
later discussion will be mostly limited to the long- wavelength regime (q — > 0). As it is 
shown in the Appendix A, the contributions from the Coulomb matrix elements with 
bond integrations [Si and/or 5 2 finite in eq. ([§])] are in this limit weak and negligible in 
comparison to the site contributions. We therefore continue by keeping only site matrix 
elements in the system (|7]), and skip <5-indices from now on. The dielectric matrix (p~0|) 



then reduces to [|T^, [UJ 



1 — VooooIIoo — ^oooi IIoi — VoiooIIio 

— Vooiorioo 1 — VoonIIoi — VonoIIio 

— rLoooIIoo — ViooiIIoi 1 — ViiooIIio 

— ^loionoo — Vionlloi — vmollio 1 



(11) 
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The discussion of the dielectric matrix in previous works [19|, |Tl|] was restricted to the 
cases in which the off-diagonal matrix elements were negligible. In the present work, we 
show that in the systems with finite interband dipolar transitions these matrix elements 
are essential and can by no means be treated perturbatively. On the other hand, the 
dielectric matrix discussed in Ref.[]19| allows for a Fermi surface crossing both bands 
due to their finite overlap. In our case of a finite gap between two bands, the matrix 
elements V" li2l ^(q, uj) on the left-hand sides of (0) are simply multiplied by di^S^i , due 
to n n (q,w) = 0. The problem is therefore reduced to a system of three equations, as is 
seen from the form of the matrix (|ilf). 

A furher simplification in [e] takes place after assuming that the products of orbitals 
<£>*(r)<£y(r) are real for each pair (1,1), as is usually true in standard cases. Then the 
matrix elements of the screened interaction (||) obey symmetry relations 

^ijMiij ^h^^h ^i-j^'i^ ' (^^) 

which are of course also fulfilled for the corresponding bare matrix elements. Note also 
that 

is generally valid due to V(r,r',u) = V(r', r, uj) and V(r + R, r' + R, u) = V(r, r', lo), 
where R is any lattice vector. 

Due to the first equality in (fL2|), first three equations in the linear subsystem (|Tl~l) 



reduce to two equations for V 0l20l > and e. g. V 0l2ll > . The corresponding matrix is given 
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by 



1 — VboooIIoo 



-vbooi(n i + ni ) 



(14) 



— V1000II00 1 — Voon(noi + IIio) 
while the microscopic dielectric function reads 

e m (q,uj) = [1 - y ooo(q)n o(q,w)]{l - Voon(q) [n i(q,o;) + IIio(q,a;)]} - 

-Vbooi(q)Kooo(q)noo(q,^) [n i(q, w) + ni (q,o;)] . 



(15) 



Finally, among sixteen TB matrix elements of the screened interaction @ six are 

given by following relations 
1 



^0000 
^0001 
^0011 

V"oioi 



[Voooo + (VoooiVboio — VooooK)oii)(noi + IIio)] 



K)001) 

1 



[Voon + (^0001^0010 — Voooo^oonjlli 



00 



[Vbioi + (VbmVoooi — VoioiVoon)(rioi + IIio)] 



Vom — — [Voin — (Vioio^oooi — VionVoooo)IIoo] 



111 



Vi 



111 



{VioioIIoo [Voioi + (VoooiVom — VoioiVoon)(IIoi + IIio) 



+Vino(n i + IIio) [Voin + (VooioVoioi — VqiuVooooJIIoo]} 



(16) 
(17) 
(18) 
(19) 
(20) 

(21) 



while the remaining ten matrix elements follow from the relations (J12|) and (|T3| ) 



The dimension of the matrix ([14]) coincides with the number of intra-band and inter- 
band transitions for the present two-band model, in accordance with the general con- 
clusion about the dimension of the dielectric matrix in the tight-binding approach from 
Ref.[|T7f| and Sect. 2. Evidently, the result ( ]I4] ) simply reduces to the one- dimensional 
matrices in the particular cases of the two-band insulator (when Hqo = 0), and of the 
one-band conductor (when Hqi = IIio = 0). 
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4 Multipole expansion 

For further considerations it is appropriate to write the bare Coulomb matrix elements 
as sums of local and long-range (two-site) parts, 

v iMtW = u iMl + w iMlW> ( 22 ) 



with 



U hk& =J d rJ ^('KCOtt^M^MO (23) 



and 



hhl , , (q) = £ e^ R fdvf drV^r - R)<(r')-^— ^; (r - R)^(r')- (24) 

Integrations in the sums (]24)) are localized within the ranges of TB-orbitals. It was 
already noted that they are small not only with respect to the range of interaction, but 
also in comparison with the distances between neighboring crystal sites. It is therefore 
natural to perform the multipole expansions of the long-range parts. This is particularly 
useful in the long wavelength limit in which the hierarchy of the leading multipole terms 
is well defined. In this Section we calculate dominant contributions to the long-range 



W, 



bare matrix elements fl24J), and indicate for which matrix elements (|22f) the on-site 
contributions (p3|) are relevant for the dielectric response in the long wavelength limit. 
Let us start from the expansion of the bare Coulomb interaction up to the second 



order in the atomic coordinates 22 



1 _ 1 R(p-p') + pp'-3(pR )(p'R ) 



p + R - p | R R 3 R 3 



1 V~V ' ' \ R 2 3i ,3 ~ 3XiXj 

o iJPiPi + PiPj) ^ " + (25) 
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where Ro = R/-R, and p and p are local electron positions at the R-th site and at the 
origin respectively. X, and pi are i-th Cartesian components of the vectors R and p 
respectively. In the next step we chose the parity of the TB orbitals by imposing 

^o(-r) = <p (r), (26) 

and 

^(-r) = ±^(r). (27) 

Since the two possibilities defined by the upper [model A] and lower [model B] sign 
in eq.(|27|) lead to qualitatively different properties of the dielectric response, we shall 
treat both in parallel. Furthermore, we choose the cubic lattice, the simplest one for 
the calculation of lattice sums in the matrix elements of the terms in the expansion 
( |25| ) . To fix ideas, we also specify in the further considerations that ipo is an s - orbital, 
while ifi is an d (upper sign) or p (lower sign) orbital, i.e. that the models A and B 
are characterized by (s,d) and (s,p) bands respectively. With these choices and with 
the assumption that there is no orbital degeneracy, the crystal symmetry is strictly 
speaking lower than cubic. Still this is not a serious inconsistency. On the one side, one 
may suppose that the Bravais lattice is built from some elongated molecules with e.g. 
a and 7r electronic orbitals and that the differences in the lattice constants along three 
orthogonal crystal directions are still small and have negligible effects on the lattice sums. 
Alternatively, the present analysis can be straightforwardly completed by including cubic 
orbital degeneracy into electronic spectrum. The dielectric matrix for the cubic crystal 
with e. g. one (partially) full s and three empty p bands will be derived in the Appendix 
B. 
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Let us now determine the leading terms in the expansion (P5|) of the matrix elements 



(22), paying particular attention to those present in the dielectric matrix ([141) . 

Voooo(q) belongs to the set of matrix elements Vi 1 i 2 i 1 i 2 (q) for scatterings in which 
electrons do not change bands. The leading long range contribution W / z 1 « 2 « 1 i 2 (q) comes 
from the monopole-monopole term which is same for all l\ and I2, no matter whether 
one (li = I2) or two (li 7^ I2) bands are involved. Since this term diverges in the 



limit q — > 0, the on-site contribution Ui x i 2 i^i 2 to the total bare interaction (|22| ) can be 
neglected. Keeping only the most divergent part of the monopole- monopole term, one 
gets for both models A and B 

47re 2 

where a 3 is the volume of the primitive cell. 

The bare Coulomb matrix elements in which one or both electrons change bands will 
depend on the choice of TB orbitals <p and ipi. The decisive difference comes from the 
fact that the matrix element for the dipole transition 

M01 = = e J <Wo(pWi(p) ( 29 ) 

vanishes in the model A and is finite in the model B. We therefore continue by considering 
each model separately. 

Let us start with the model B. The leading long range contribution to the matrix 
elements in which one electron changes the band comes from the monopole-dipole term 
in the expansion (^), and reads 

WV(q)=E^ R ^^- (30) 
R^O U 
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(I ^ I). The divergent sum in this expression depends on the angle between q and n U ', 
where the direction of /j, u ' is determined by r - dependences in ip and tp\ (|29|). With 
if i chosen to be the p - orbital elongated along the x - axis, fx 01 is also directed along 
this axis. Noting in addition that in the model B the on-site term U m > vanishes due to 
symmetry reasons, one gets 

Vuu'(<l) = -Vi'm(^) = — § 2 31 

where \x = |a* 01 |. The leading long range contribution to matrix elements with both 
electrons changing bands is the dipole-dipole one, which after performing the dipolar 



summations 23] for the cubic crystal lattice in the limit q — > 0, reads 



^Y(q) = ^(^-l)+^Y, (32) 

with I 7^ I . Eqs.(^), (|31| ) and (|32|) exhaust all matrix elements for the model B. We 
conclude that for q — > the on-site and the long range contributions enter on the equal 
footing only into the matrix elements (0), while in the matrix elements ( p8| ) and ( |3T| ) 
the former are negligible with respect to the latter. 

In the model A the leading long range contributions to the matrix elements V ull > (q) 
and V^';'(q) come from monopole-quadrupole and quadrupole-quadrupole terms in the 
expansion (|25|) , respectively. If the matrix elements of the quadrupolar transition be- 
tween the orbitals (fi and (pp are finite, these contributions behave in the limit q — > as 
q° and q 2 , respectively. We do not need their detailed forms, since the further conclu- 
sions will be based solely on the fact the corresponding total matrix elements V ml >(q) 
and Vyjy(q) are regular in the long wavelength limit (and V^y ~ £^«y i n addition). 

The most important difference between the q — > limits of models A and B appears 

1G 



in processes in which one electron changes its band. While in the model B the corre- 
sponding matrix element (31) comes from the long range part and diverges as the 



analogous matrix element in the model A is the sum of the long range part and the 
on-site contribution and behaves as q°. 

It remains to determine the asymptotic behavior of the polarization diagrams which 



figure in the dielectric matrix ([14]). We limit the present discussion to the regime (q — > 
0, uj 0) in which the real parts of the polarization diagrams are given by the well-known 
expressions 

n ^ 

Re n 00 (q, uj) = — e - — , (33) 
m* uj 2 

and 

Re[U 01 (q,uj) + U 10 (ci,u)] = 2n e Eg (34) 

uj z — hi* 

m* is the effective mass and n e is the number of electrons per site in the lower band, while 
Eg = E\ — Eq is the energy difference between the band centers. Here we neglect the 
bandwidths in the expression (|34|) assuming that they are much smaller than E g . Note 
that the expressions (^) and ( |3"4] ) do not depend on the symmetry of TB orbitals <p and 



(pi. Furthermore, we do not introduce the corresponding imaginary parts, relegating the 
discussion of the effects of Landau damping to the paper II. We only note that in the 
above (q, c<j)-regime Im LToo = while Imflloi + ITio] may be finite. 

We are now ready to compare the microscopic dielectric functions for models A and 
B. Let us start with the simpler case of model A, already considered in Ref.([l9|]) and 
explored in Refs.(p(| ^TJ). As it is seen from (|33D and (j34"D, the term in the expression 



15|) which comes from the off-diagonal elements in the matrix (|Tj) vanishes in the long 



17 



wavelength limit. In other words, the intra-band and inter-band polarization processes, 
given by the diagonal elements in this matrix, become decoupled, and it is legitimate 
to distinguish between the intra-band (ej ntTO ) and inter-band (ej nter ) dielectric functions 



19| , the total dielectric function being the product of the two. Due to this factorization 
the collective modes, defined by the real zeros of e m in the u;-plane, are also either 
of intraband or of interband origin. They will be analysed in detail in the paper II. 
Furthermore, the expressions ( |T6| - |2T|) for the screened Coulomb interaction reduce to 

Va'W /Vu'W = l/ei n tra (35) 

(with either I = I' or I ^ I'), and 

Vuw /Villi' — Viii'i'/Viu'i' = l/e inter (36) 

(with I 7^ /'). Thus all purely intraband scatterings are screened by ei ntr a, while the 
scatterings with at least one interband electron transition are screened by Winter- 

The above decoupling does not hold in the model B. The product of two off-diagonal 
matrix elements in (JT^) now remains finite in the long wavelength limit, and introduces 
a coupling between the intraband and interband contributions into the microscopic di- 
electric function (|I^). The present analysis shows that this coupling comes from the 
divergent monopole-dipole contribution ([IT]) to the matrix element V^/(q). Its most im- 
portant consequence, the hybridization of the intraband and interband collective modes 
in conductors, will be considered in the paper II. 
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5 Macroscopic dielectric function 

The aim of this section is to establish the relation between the macroscopic dielectric 
function ejy(q, uj) and our microscopic dielectric function e m (q, uj). The former is defined 
in a standard way as the ratio between the Fourier components of the spatial averages 
of the bare (external) and screened (total) Coulomb potential. To get macroscopic 
averages, it is enough to take the unit cell as the averaging volume 0. Let us denote 
the averaged bare and screened interactions between the unit cells at sites R and R' 
by V av (H, R') and V av (H,H',uj) respectively. Since at long wavelengths details in the 
space dependence of the true (non-averaged) interactions on the scale of unit cell are 
irrelevant, the true and averaged interactions are in this limit indistinguishable, and one 
gets the connection 

V av (q + G, q + G', uj) = S G , 5 G r ;0 V(q, q, uj) (37) 

between the corresponding Fourier transforms. Here G and G' denote the vectors of the 
reciprocal lattice. The same relation holds for the averaged bare interaction V av (q). 
The macroscopic dielectric function is now defined by 

— 7 \ = — TTT~\ — = e fa> u • 38 

e M (q,u) V av {q) 

The ratio V av (q, q, uj)/V av (q) can be also interpreted as the q-component of the total 
charge density, redistributed due to the Coulomb screening, averaged across the unit 



cell and divided by the corresponding q-component of the probe charge [|17]]. The last 
equality in eq.(^) follows from the standard definition of the dielectric matrix in the 
plane wave representation, V(q, G, G', uj) = e _1 (q + G', q + G, u>)V(q, G, G), and the 
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relation (|37l) . Our aim is however to express ej\,/(q, uj) in terms of the TB bare and 
screened Coulomb matrix elements (|22|, |8|). The TB matrix elements (§) can be expressed 
as double sums over the Fourier transforms V{<\ + G, q + G', u) through the relation 

v hl2VA ( q , W ) = £ v (q + G ) J ^ 2 (q + G Wq + g, q + g', w ) (39) 

G,G' 

with 

I«,(q) = |rfr^(r)e^(r). (40) 
The corresponding relation for the averaged matrix elements follows after inserting the 
relation ([37]) into fl39|) . Expanding the matrix elements fl40|) in powers of q we get in the 
limit q — > 



zqu, ,/ 



y a «(q,^) (41) 



e 

with /a^j/ being the dipolar matrix element (p9|). Thus, whenever ^ 7^ i x and/or l 2 ^l 2 , 
the corresponding ratios ^^^^'(q, a>)/V ,;(q, a;) tend to zero as q — ► 0. In other 
words, since the averaged interactions do not vary within a unit cell, the contribu- 
tions from the dipolar and higher transitions to the corresponding TB matrix elements 
V av i^ff (q, oj) vanish. For l\ = l[ and l 2 = l 2 , i-e. for matrix elements with a finite 
monopole-monopole contribution, we have 

V av , h i 2h i 2 (q,u) = V hhh i 2 (q,uj) = V av (q,u) (42) 

for any pair of orbital indices l\ and l 2 . The first equality in ( f4"2| ) is based on the same 
argument as the relation (PT|). Since the relations ( |4"T| ) and ( ^2]) are also valid for the 
bare interaction V(r — r'), the macroscopic dielectric function reads 

e M (q,w) = = - ( r- (43) 
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This result is generally valid for any number of bands. It also does not depend on 
the method of calculation of the linear dielectric response, which brings in a particular 
relationship between the screened and bare Coulomb interaction. The independence 
of the right-hand side of fl43f) on the indices li and Z 2 may serve as a check of the 
consistency of a given approximation in the macroscopic limit. It is easy to see that this 
independence is realized for the RPA result ( |T6| - |2T| ). 



Although e^(q, to) is the ratio of the bare and screened monopole-monopole inter- 
actions, this of course does not mean that other higher order terms in the multipole 
expansion do not contribute to the macroscopic response. This becomes clear already 
for the simple two-band models from Sects. 3, 4. For the model A the result (|35| ) reads 



ejif (q, w) = e intra (ci, u) = 1 + 4vra c (q, u), (44) 

with 



e 2 



« c (q, w) = -^^n 00 (q, uj) (45) 

representing the intraband polarizability. Thus, the macroscopic dielectric function is 
entirely determined by the intraband processes, and reduces to unity when the lower 
band is fully occupied. In other words, the interband modes are not directly optically 
active, and can possibly be observed only by Raman scattering. 

The macroscopic dielectric response for the model B follows after inserting any of the 
expressions (|I^JT9|j2l|) into (|43D , and taking into account the long wavelength asymptotic 



expressions (|28|j3l| , |32|) . One gets 

Cm _ 1 47T/i 2 

6m 3a 3 



(n 01 + n 10 ). (46) 
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Generally, the right-hand side in this expression depends on the direction of q due to 
the anisotropy of the polarization diagram (Hoi +II10). However, in the particular limit 
given by eq.(|34"D this is not the case, so that the orientational dependence of 6Af(q, uj) 
is determined solely by that of the microscopic function e m (q, uj). Let us consider here 
only two particular directions, q _L /x and q || /x. 

For q = q^ _L /x the bare monopole-dipole matrix element fl3T|) vanishes, so that 
e m = tintratinter with e in t er equal to the right-hand side od (pEE|). Hence one gets 

£A/(q±, w) = l + 47ra c (q_L, uj), (47) 

as should be expected, since for q _L \x the interband longitudinal polarization processes 
are not possible. 

For q || /x the (intraband) monopole - (interband) dipole coupling is finite and, in 
addition, even for an insulator there is no simplification in (f46|) as for q _L /x. The 
macroscopic dielectric function can be then written in the form 

e M {q h oj) = 1 + 4ttq! c + - ^ =z . (48) 

1 3~^-f 



Here 



a Mb u ) = TTTt 37 V7^' 49 



with ai(q\\,u) being the dipolar interband polarizability 



a/(?||,w) = — - n i(g||,a;) +U 10 (q h uj)}. (50) 

(X 



The three quantities which enter into the expressions (f|^, [50D represent the intraband 
(a c ), interband («/) and intra-atomic [C/oon(noi + n 10 )] contributions to the macro- 
scopic dielectric response. Consequently, this expression covers, and interpolates be- 
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tween, limiting cases of one-band conductor, two (multi)-band insulator, and atomic 
(zero bandwidth) insulator. Before going into particular limits, we remind that Adler 
parametrized any deviation of the expression (f|8j) from the usual form containing 



only the dipolar interband polarizability floTf ) by introducing a self-polarization correc- 
tion, defined as aj(q\\,u) — aj(q\\,u)) = C. When expressed in terms of TB quantities, 
this parameter reads 

c _n 2 twrioi + n 10 ) 2 

As it will be shown below, the TB results (^) and flSDp provide a direct physical insight 
into this correction for atomic insulators as well as for band insulators and conductors. 
Let us now consider some characteristic limits. The usual Lindhard (i.e. Sellmeyer 



24] ) expression for a metal follows after neglecting interband polarization (aj = 0), while 



the opposite limit for an atomic insulator is obtained after putting a c = and neglecting 



all bandwidths The latter result is the well known Lorentz-Lorenz (i.e. Clausius- 
Mossotti) expression. In this limit the polarizability aj which enters into eq.(^8[) is given 
by 

aiM = - — ^ _ E2 J^ , (52) 



i.e. it does not differ in form from the expression ( |5T1|J5^ ). Here [i^ff = V Eg/Eg,eff an d 
Eg e ff = Eg + 2nE g UoQu are the renormalized values of the dipole matrix element and 
the level spacing, i.e. instead of the initial (e.g. atomic Hartree-Fock) parameters one 
obtains those which include the on-site RPA screening. In other words, the Adler's self- 
polarization correction in this case simply reduces to a renormalization of the parameters 
in the dipolar (intra-atomic) polarizability. This reflects the fact that the RPA scheme 
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proposed in the present work treats both nonlocal and on-site screenings at the same 
level of approximation. In this respect, we warn that in the limit (|5^) the RPA is usually 



not a reliable scheme at the atomic (or molecular) level [25], and that it is obviously 
invalid when the on-site interactions are large in comparison with E g . 

In the cases of insulators and conductors with finite bandwidths the effective dipolar 
polarizability (f£9|) does not have the simple form and, moreover, cannot be reduced 



to the form of an effective interband polarization diagram, cf. (|50|). Instead, one has a 
nontrivial a;-dependence in the denominator of the expression (pEDp. The self-polarization 
correction is no more of entirely local (i.e. intra-atomic) origin, but appears to be a 
combined effect of the on-site Coulomb interaction L/0011 in which two electrons exchange 
their orbital states, and of the finite inter-site tunneling (i.e. of the finite bandwidths). 



The role of £/oon in the macroscopic dielectric function was already noticed in Ref.jl l 



We note that some later results for eu 0; Hi) although based on the RPA method 



cannot be expressed in the form (|49| , p0[ ) i.e. ([51]). The reason might be traced in the 
additional approximations which, unlike the present RPA scheme, do not include the 
local on-site screening on the same footing with the long range one. 

Furthermore, it is appropriate to point out that aj appears in the denominator of 
eq.([ISp due to the finiteness of G 7^ terms (local fields) in the Fourier transform of the 
dipole-dipole interaction d32|) . The G^O terms are comparable to the G = contribu- 
tion and should not be neglected, as already stressed in the Introduction. Their omission 
transforms the formula fl48f ) into an expression which is additive in the polarizabilities 
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6 Conclusion 

The most important step in the present work is the use of the complete set of TB 
states in the representation of the bare and screened Coulomb interactions. In this 
basis, the symmetry properties of the Coulomb matrix elements follow directly from 
the symmetry of molecular orbitals participating in the band states. As a consequence, 
the initial system of linear equations for the screened matrix elements decomposes into 
smaller subsets with a common dielectric matrix. E.g., in the two-band case sixteen 
linear equations were reduced to two equations. This agrees with our general conclusion 



171 that for a multi-band system with only two-center Coulomb interactions taken into 
account, the TB dielectric matrix within RPA has a dimension equal to the number of 
non- vanishing intra-band and inter-band polarization diagrams (H). 

The present approach is particularly convenient in the long wavelength limit, in which 
the multipole expansion is feasible. Again, the symmetry of molecular orbitals is of 
central importance, since it determines the degree of mixing between intra-band and 
inter-band polarization processes in multiband conductors. In particular, whenever the 
interband processes are dipolar the off-diagonal elements in the TB dielectric matrix are 
finite and cannot be treated perturbatively. The inclusion of the long-range monopole- 
dipole interaction induced by these processes leads to the correct expression for the 
macroscopic dielectric function, and, as will be shown in paper II essentially influ- 
ences the spectrum of collective modes for a multiband electron liquid (see also Ref. 



ITU). 111 this respect it is important to recognize that the results for e m (q, u) ( |10D and 
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ejvf(q, oj) (fP|) reduce to the additive form || H 

e(q, to) — 1 = 47ra; c + 47T a/, (53) 

r 

only in the high frequency limit, i.e. for frequencies much higher than all interband 
energy differences E g i. [The /-summation in eq.(|53"D goes over all relevant interband 
polarizabilities ([SOD]. It is interesting to note that in this limit we apparently recover 
the well-known expression for the high frequency plasma edge expressed through the 
standard plasmon sum rule, involving the free electron mass. 

In conclusion, we point out that the present approach might have numerous applica- 
tions, particularly in systems with dipolar degrees of freedom which are sensitive to the 
variations of the valence band properties induced by pressure, changes in band filling, 
etc. We expect that the present method will prove very useful in such investigations. 

A Appendix A 

In this Appendix we consider the influence of the site-bond (<5i or 82 different from zero) 
and bond-bond (Si and 82 different from zero) Coulomb interaction (|8J) on the dielectric 
response in the long-wavelength limit (q — > 0, uj finite). We take the simplest examples 
of one-band and two-band systems with the one-dimensional lattice [Si in eqs. (|7| - |T0D 
equal to and ±1]. The extensions to more complex band structures and lattices are 
straightforward. 

The system (|7]) simplifies after assuming, like in Sects. 3 and 4, that the orbitals ipi(r) 
are even or odd, and that the products ip* (t — S)ipi 2 (r) are real. This leads to additional 
symmetry relations between bare Coulomb matrix elements (|8|), which in the case of 
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two-band system = 0, 1) read 

V hhl i A (5 u 6 2 ) = (-1)£^ 2 ^(<5 2 A) = (-i)£k Vm2h{ _ 5 ^_ 5l) = 

= Vw^idt, -6 2 ) = Vi'^i-St, 6 2 ) = V mih (-8 X , -S 2 )- (A.l) 

Due to these symmetry relations, the number of equations in (|7[) reduces to N e /2 or 
(N e + l)/2 for even or odd N e respectively, where N e is the product of Z + 1 and the 
number of the non-zero polarization diagrams among n o, n i and n 10 . 

In particular, taking into account these relations, we obtain the microscopic dielectric 
function for one-band metal in the form of a 2 x 2 determinant 

1 - E$ V{0, S)U (S) - Z$ V(l, S)U Q {S) 

- E S V{0, S) [U (S - 1) + n (<5 + 1)] 1 - Es V(l, S) [U (S - 1) + U (d + 1)] 

(A.2) 

with V(<5i, 62) = Voooo(^i, ^2, q) and n (<5) = n o(<5, q, oj). In the limit q — > the 
leading monopole-monopole terms in the bare matrix elements appearing in ( |A.2j ) are 
simply related by 

V(l, 1, q) = SV(1, 0, q) = S 2 V(0, 0, q) = S 2 — , (A.3) 

a 6 q z 

where S is the direct overlap After inserting these relations and the expressions (0) 
for the polarizability diagrams n (0),IIo(±l) and n (±2) into the expression (|A.2p , it 
reduces to 

e m = l- V{0, 0, q)n (0) [l + 0{q 2 )] . (A.4) 

Thus, we conclude that in the limit q — >• 0, u finite, the microscopic dielectric function 
for a one-band metal does not depend on the overlap factor S which enters into the 
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site-bond and bond-bond interactions through the expressions ( |A.3|) . In other words, it 
is completely determined by the site-site correlations. 

In the case of the two-band insulator the most important simplification comes from 
the fact that the only non-vanishing polarization diagrams are the inter-band ones with 
the zero phase factor [8 — d' = in eq.(S)]. The microscopic dielectric function then 
reduces to 

1-^(1,1)11! -^(0,1)11! -^(-1,1)11! 

-17(0,1)11! 1-^(0,0)11! -^(0,1)111 (A.5) 

-y;(-i,i)n! -Vi*(o,i)ni 1-^(1,1)^ 

with Vi(S,S') = Viooi(6 , 6' ; q) . Since only one polarization factor [n i(0,q, a;) + 
IIio(0, q, uj) = III] figures in this determinant, the criterion of weakness of site-bond 
and bond-bond correlations with respect to the site-site ones obviously involves only 
the corresponding matrix elements of bare Coulomb interaction, i.e. Vi(0, 1), Vi(l,l) 
and Vi(— 1, 1) vs Vi(0, 0). More precisely, the three- and four-center correlations may be 
considered as weak perturbations provided the matrix element of bond dipolar transition 



fi(S) = J dr^*(r)rv?i(r - 6) 



(A.6) 



and local site-bond and bond-bond matrix elements are much smaller than the respective 
site parameter fl29[) and t/ooii [see eq. (|32|) l. 

For a two-band metal e m is given by a 5x5 determinant, not written here for the 
sake of space. Although this determinant contains intraband and interband polarization 
diagrams (^) with various phase factors, the analysis of the long wavelength limit leads 
to the same conclusion, namely that the criterion for the weakness of site-bond and 
bond-bond correlations in the dielectric function involves only the corresponding bare 
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Coulomb matrix elements. 



B Appendix B 

In Sections 4 and 5 we have assumed within the model B that there is one orbital with 
p x symmetry at each site of a cubic Bravais lattice. In order to complete the discussion, 
we take here the simplest generalization by which the electron bands have the cubic 
symmetry. We assume that there are one lower s-orbital ( |2"6D and three degenerate 
upper p - orbitals with the symmetry 

<P*(-x, V, z) = -<p x (x, y, z) (B.l) 

and equivalent ly for x — > y — > z. 

Following the procedure from Sec. 3 we get for this four-band model a linear Dyson 
system (0) which contains 256 equations and is decomposed into 16 subsystems of 16 
equations. After assuming, like in Sec. 3, that all products <^*(r)<^/(r) are real, the 
problem reduces to a system of four equations 



^OOOx 


— VqoOx + K)000 ^000x1^00 


+ K)00x^00a;a;(n0i 


+ n x0 ), 


(B.2) 




= Vofte + KrOOO^OOOxnoo 


+ VoOxxVoOxx(Jl 0x 




(B.3) 


^OOj/x 


= VoQyyV 00y X (Tl0y + Hyo) 






(B.4) 










(B.5) 



Here subscripts x, y and z stand for p x ,p y and p z orbitals respectively, and the wave 
vector q is oriented along the crystal axis x. Due to the isotropy of the cubic crystal this 
special choice does not harm the generality of the system ( |B.2[ - |B.5| ) (see also Ref. |T] 
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The corresponding dielectric matrix is given by 



1 — V0000IT00 — VoOOa;(no2; + ^-xo) 

— V^oooIIoo 1 — Vboxx(no»c + II x o) 

o o l-VoowCnoy + nii/o) 

o o o i-v 00zz {u 0z + u z0 ) 



(B.6) 

and all TB matrix elements V^^y with = 0,x,y,z follow straightforwardly, but we 
do not write them in order to save space. 

The dielectric properties of the isotropic crystal in the direction parallel and per- 
pendicular to the the wave vector are described by the upper left 2x2 block, and by 
the third and fourth diagonal element in the matrix (|B.6|) , respectively. Note that the 



parallel and perpendicular responses with respect to any direction of the wave vector 
in the cubic crystal are identical to those along the respective preferred dipolar axis x, 
and the perpendicular axes y and z within the model B. Thus, all conclusions derived 
within the model B for the preferred direction x are valid also for any direction in the 
isotropic cubic crystal. In particular the expression (^) for £j^(q, o>) is now valid for 
any orientation of q. 
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